Landau Cooling in Metal-Semiconductor Nanostructures 
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An electron-cooling principle based on Landau quantization is proposed for nanoscale conductor 
systems. Operation relies on energy-selective electron tunneling into a two-dimensional electron gas 
in quantizing magnetic fields. This quantum refrigerator provides significant cooling power (~ 1CP 9 
W at a few K for realistic parameters) and offers a unique flexibility thanks to its tunability via the 
magnetic-field intensity. The available performance is only marginally affected by nonidealities such 
as disorder or imperfections in the semiconductor. Methods for the implementation of this system 
and its characterization are discussed. 
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The relevance of understanding and tailoring thermal 
transport in nanoscale and mesoscopic conductors is be- 
coming more and more apparent P, ^ S 0] • This stems 
from the ever increasing miniaturization of the solid- 
state systems of interest for optoelectronics and advanced 
condensed-matter research which in turn demands a pro- 
found understanding of the generation and manipulation 
of heat flow. In this context electron cooling is prob- 
ably one of the most relevant aspects. Any solid-state 
system characterized by a gap in the density of states 
(DOS) at the Fermi energy is, in principle, suitable for 
refrigeration schemes In fact a gap provides an 

energy-selective electron-filtering mechanism which only 
allows hotter quasiparticles to be extracted efficiently. 
Superconductors represent prototype systems for refrig- 
eration schemes operating around or below IK [H,0|, and 
have yielded substantial electron temperature reduction 
in metallic nanostructures [f|. 

Here we propose a novel mechanism for refrigeration 
based on Landau-level quantization. To this end we ad- 
dress thermal transport in two-dimensional electron gases 
(2DEGs) showing that heat can be efficiently extracted 
from a metallic nanostructure, thanks to the peculiar bulk 
DOS originating from Landau quantization. We find a 
large heat current through the system (~ 10 -9 W) yield- 
ing a sizable electron cooling also at bath temperatures 
as high as a few K. Moreover, the possibility to easily set 
the magnetic field value makes it possible to finely tune 
and optimize the performance available with this refrig- 
eration principle under different operational conditions. 

The system we envision is sketched in Fig. 1(a) and 
consists of a normal-metal island (N) coupled to the bulk 
of a wide 2DEG region through a tunnel barrier of re- 
sistance i?t, in an open pseudo Corbino-disk geometry 
0, [H • A magnetic field B is applied perpendicularly to 
the plane of the 2DEG, and the whole system is biased 
at voltage Vb through a superconducting lead (S) in good 
metallic contact with N to prevent heat injection [l[ . In a 
quantizing magnetic field the bulk DOS T>{E) of a 2DEG 
in the presence of impurities or disorder can be expressed 
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where E is the energy measured from the Fermi energy 
E F , C = 1/(2tt2V|), V = m*/{2irh 2 ) is the DOS per 
spin for B = 0, £b = \J h/\eB\ denotes the magnetic 
length, r = y / 2h 2 uj c / (tit), and E° = (n+^)huj c -^af>ioz- 
Furthermore, lo c = eB/m* is the cyclotron frequency, 
fiwz = g*/iBB is the Zeeman splitting, m* and g* are 
the electron effective mass and the effective g-i actor, re- 
spectively, r is the relaxation time for B = 0, /ib is the 
Bohr magneton, and a = ±1 is the spin index. As we 
shall demonstrate below, the existence of Landau gaps 
will lead to a significant heat flow through the system 
thus implying Zioi-excitation removal (i.e., refrigeration) 
from the N island. Within a phenomenological tunnel- 
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FIG. 1: (color) (a) A normal metal island (N) is tunnel cou- 
pled to a two-dimensional electron gas (2DEG) in an open 
pseudo Corbino disk geometry in the presence of a quantizing 
magnetic field. Gray region represents an ohmic contact to 
the 2DEG. (b) Heat current Q extracted from the N region vs 
voltage drop V across the barrier, calculated for different fill- 
ing factors at k^T = 10 _2 -Ef. Other parameters used in the 
calculation 10]: TE ¥ /(2ixh) = 484 and lu z /lu c = 14.7 x 1(T 3 . 
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ing Hamiltonian approach, the heat current Q(V) flowing 
out the N side of the 2DEG/N junction can be written 
as Q 



Q(v) 



1 
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dEEV(E)[h(E) - / 2D eg(£)], (2) 



where E = E— eV. V is the voltage drop across the bar- 
rier, and /n(2DEG) is the Fermi distribution at tempera- 
ture T in the N island(2DEG). Furthermore, V = V -V th 
where Vth is a threshold voltage 14]- The heat current 
Eq. (|2|) can be rewritten as a series expansion in powers 
ofr/(fc B T): 
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where A° m = ^n(E% - Ep - eV) if m is even, A" m = 

{m -l)\\ l if 171 is ° dd ' and /n(2DEG) = 9 "7n(2DEG) / 9E m . 

Equation ([3]) simplifies considerably the calculation of 
QiV) in the high-temperature limit, where only a few 
terms need to be retained. 

Figure 1(b) shows the heat current Q(V) versus bias 
voltage V at k B T = 10^ 2 Ep calculated for several v val- 
ues [lOj. Q is a symmetric function of V. For even v it 
presents a positive peak (describing heat extraction from 
the N region) in correspondence of the Landau level clos- 
est to the Fermi energy, and thereafter rapidly decreases 
to negative values. Furthermore, the maximum value of 
Q decreases by increasing the filling factor, reflecting the 
DOS-amplitude reduction at smaller magnetic fields. For 
odd values of v [see the case of v = 1 in Fig. 1(b)] the 
heat current is always negative, unless Zeeman coupling 
completely lifts spin degeneracy giving rise to a gap [ll| . 
This full tunability of the heat current is a very attrac- 
tive feature readily accessible by setting the magnitude 
of the applied magnetic field. The impact of disorder in 
the 2DEG on heat transport is analyzed in Fig. 2 (a). 
Q(V) at k B T = 10~ 2 E F and v = 4 is calculated for sev- 
eral values of the relaxation time r. The latter is related 
to the 2DEG mobility, /i, according to /i = er/m*, and 
describes (with the chosen parameters) a realistic GaAs 
2DEG with mobility values in the range between ~ 27 
and 1300 m 2 /Vs. In particular, by shortening the re- 
laxation time (i.e., by increasing disorder) the positive 
Q peak broadens and its amplitude is suppressed with 
the maximum moving towards lower bias voltage. This 
stems from disorder-induced broadening of the DOS [see 
Eq. (UJ)], and indicates that high- mobility 2DEGs are 
required to maximize Q. We stress, however, that a re- 
duction of t by a factor of 5 (e.g., from t = 5 x 10~ 10 s 
to 1 x 10~ 10 s for a GaAs 2DEG) yields only a small sup- 
pression of Q (~ 14%), and that long relaxation times in 
2DEGs can be achieved nowadays with epitaxial growth 
techniques: relaxation times as long as 2 x 10~ 10 s (which 
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FIG. 2: (color) (a) Heat current Q vs V calculated for several 
scattering times r at v = 4 and feT = 10 _2 i?F. (b) Q vs V at 
v = 4 calculated for different temperatures and tEf / (2irh) = 
484. The Zeeman splitting is luz/uJc = 14.7 x 10" 3 . 



corresponds to /i ~ 500 m 2 /Vs) and even larger can be 
obtained in the GaAs/AlGaAs material system 12, 13j |. 

The role of temperature on thermal transport is dis- 
played in panel (b) of Fig. 2 where Q vs V at v = 4 
is shown for several values of T. At each tempera- 
ture Q is maximized at an optimal bias voltage (Vopt) 
which, at the lowest temperatures, approaches the sepa- 
ration between the Fermi energy and the closest Landau 
level. By increasing T, the heat-current amplitude grows 
and reaches its maximum Q ~ 5 x 10 _3 £'|/(e 2 i?t) at 
V^pt ~ 0.185£ F /e for k B T = 5 x lO^p. A similar 
qualitative behavior also occurs at the other filling fac- 
tors. At higher temperatures Q becomes negative. The 
behavior of the maximum heat current Q max (i-e., that 
at the optimal bias voltage) as a function of temperature 
is shown in Fig. 3(a) for different filling factors v. In par- 
ticular, at low T, Qmax is a bell-shaped function which, 
depending on v, is maximized at different temperatures. 
The optimal temperature decreases by increasing the fill- 
ing factor while, at the same time, the window of positive 
heat current is reduced. For instance, at v = 2 the heat 
current is maximized at k B T ~ 0.1 Ep where it reaches 
the value Q^= 2 ~ 2.1 x 10~ 2 £|/(e 2 i?t) while, at v = 4, 
Q v =^ ~ 5.3 x 10- 3 £|/(e 2 i?t) at k B T ~ 0.045£: F . By 
choosing R t = 500 Q the maximum heat current corre- 
sponds to Q^= 2 ~ 4 nW and Q^_^ 1 nW at T ~ 11.6 
K and T ~ 5.2 K, respectively [lfj. Furthermore, for 
k B T > 0.2Ep, Q maK (T) increases. This stems from the 
contribution to the heat current brought by additional 
Landau levels that occurs at higher temperatures. 

Whether a cooling principle is usable in practice de- 
pends, in general, on its thermodynamic efficiency (rj), 
defined as the ratio of heat extracted from the system 
and the work done by the source, rj{V) = Q(V)/[VI(V)] 
0, where I(V) = ^ dEV(E)[f N (E) - / 2 deg(£)] 
is the charge current flowing through the 2DEG/N junc- 



tion. The function 77 versus T calculated at V op t is 
shown in Fig. 3(b) for the same filling factors of Fig. 
3(a). At low temperatures, r\ is a non-monotonic func- 
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FIG. 3: (color) (a) Q max calculated at V op t vs T for various 
filling factors. In the low temperature regime the heat current 
is maximized at different temperatures depending on v. (b) 
Efficiency 77 vs T at the same v as in (a). Other parameters 
used in the calculation: tEf / '(2nh) = 484 and wz/w c = 14. 7x 
10" 3 . 
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tion of T and presents a relative maximum taking val- 
ues as high as ~ 27% almost independently of v (note 
that such a value for rj is comparable to the efficiency of 
a static superconductor-insulator-normal metal electron 
cooler which is of the order of ~ 25% [l|). After that, 77 
starts to increase monotonically at a temperature which 
depends on v, and saturates asymptotically around 60% 
for large T. As we shall show in the following, however, 
only the lower-temperature behavior is relevant for elec- 
tronic cooling. 

We now study the consequence of the positive heat cur- 
rent flowing off the normal metal, i.e., electronic refriger- 
ation of the N island. We focus on the quasi-equilibrium 
limit [l[ which can be attained at sufficiently low tem- 
peratures (i.e., for metallic systems typically around or 
below 1 K) when fast inelastic electron-electron relax- 
ation drives the quasiparticle distribution toward a Fermi 
function at an effective temperature (T ) which can dif- 
fer significantly from that of the lattice (Tbath)- This re- 
quires the tunneling rate of electrons in N to be (much) 
smaller than the electron-electron relaxation rate (t~*). 
We note that in pure metals the latter is of the order 
of t~1 ~ 10 8 . . . 10 10 s _1 [l]| and even faster, so that 
the quasi-equilibrium regime turns out to be a realistic 
assumption for our system In order to investigate 
the electronic refrigeration of the N island, we need to 
evaluate those mechanisms which tend to drive power 
into the N electrons. The main heat load comes from 
electron-phonon interaction with energy exchange given 
by [3 Qo-bath = SV(T 5 -T b 5 ath ), where E is the material- 
dependent electron-phonon coupling constant and V is 
the volume of the N region. Here, we neglect further 
sources of heat due to other relaxation mechanisms, for 
instance, radiative electron-photon heat load [3, [13], an d 
we suppose negligible thermal-boundary (Kapitza) re- 
sistance between the N phonons and the 2DEG lattice 
[E EH ■ The stationary T G reached by the N island elec- 



FIG. 4: (color) (a) Temperature T e of the N island vs V 
for different filling factors at Tbath = 1-5 K. (b) T c vs V for 
v — 4 at a few bath temperatures, (c) Minimum temperature 
To, min of the N island vs v at a few bath temperatures. The 
curves correspond to a realistic GaAs 2DEG with Ef — 10 
meV, r = 2 x 10~ 10 s, uz/w c = 14.7 x 10~ 3 , and an aluminum 
island of volume V = 5 x 1CT 21 m 3 , E = 2 x 10 s WK^m" 3 , 
and Rt = 500 Q. 

trons is obtained by solving 

Q(V, T e , Tbath) + Qc -bath 

(T , Tbath) =0, (4) 

which requires the sum of all thermal fluxes in the N 
island to be zero. In writing Eq. (01 we set the elec- 
tron temperature of the 2DEG equal to Tbath [H, which 
is justified by our assumption of a macroscopic 2DEG 
reservoir where electron-phonon interaction yields com- 
plete thermalization. 

In Fig. 4(a) we plot the steady-state electron tempera- 
ture T e in N versus bias voltage V at lattice temperature 
Tbath = 1-5 K for a few filling factors. The electronic 
temperature is obtained as a solution of Eq. assum- 
ing an aluminum (Al) islan d 1 19j of volume V = 5 x 10~ 21 
m 3 , S = 2xl0 8 Wr 5 m- 3 [Jand R t = 500 Q. For each 
v, the curves demonstrate a drastic reduction of T e with 
a minimum localized around the corresponding Landau 
level. Furthermore, the minimum of T e depends on v, i.e. 
on which Landau level is more effective at fixed T bath in 
removing hot excitations from the N island. In partic- 
ular, at v = 6, a minimum T e of about 228 mK can be 
achieved, which corresponds to temperature reductions 
of the order of 85%. We emphasize the magnitude of 
this value despite the large electron-phonon interaction 
present at such high lattice temperatures. Figure 4(b) 
displays T c versus V at v — 4 for different bath temper- 
atures [T e (V = 0) = Tbath]- The results show that a sub- 
stantial electron temperature reduction can be achieved 
over a broad range of Tbath- At the lowest Tbath the 
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temperature reduction becomes less pronounced. Never- 
theless, the magnetic field can be tuned in order to opti- 
mize the cooling process even at the lowest bath temper- 
atures. The minimum achievable electron temperature 
(7e,min) versus v (or B) is shown in Fig. 4(c) for three 
different bath temperatures. T^mm shows an oscillatory 
behavior as a function of v which can be understood as 
follows. For odd u, no cooling is possible owing to the ab- 
sence of a Landau gap at the Fermi energy. Consequently 
T^min = Tbath giving rise to a maximum. r Cimin then 
takes its minimum values between two consecutive odd 
values of v. We also note the overall negative slope of the 
7e,min(^) curves in the regions between consecutive odd 
v values. This fact indicates a reduction of the minimum 
achievable temperature at low B, which originates from 
the larger cooling efficiency provided by smaller Landau 
gaps at these bath temperatures [see Fig. 3(b)]. 

We should like to further comment on the experimental 
realization of this system. A 2DEG in the GaAs/AlGaAs 
material system is a suitable candidate to exhibit Lan- 
dau quantization although semiconductors with 
larger (/-factor (such as InAs 2DEGs) could provide ad- 
ditional flexibility thanks to enhanced magnetic-field- 
induced Zeeman splitting. Furthermore, Fermi energy in 
the 2DEG can be changed through electrostatic gating, 
thus providing an additional degree of freedom for tuning 
the electronic temperature at different magnetic fields. 
As far as the S biasing contact is concerned [see Fig. 
1(a)], a NbN electrode is particularly suitable to inject 
electric current in the N region. Its upper critical field 
can exceed 20 T at routinely accessible temperatures, 
thus allowing one to reach low v at high magnetic fields 
without suppressing superconductivity [2(|. In practice 
T c can be probed in the N island by Coulomb blockade 
thermometry [211 ] which proved to be virtually immune to 
high magnetic fields [13] ■ A shot-noise thermometer [23| 
is a further possibility. Finally, we point out that the 
present system is to a large extent unaffected by non- 
idealities such as barrier nonuniformity which typically 
degrades cooling with a superconducting gap [l|. 

In summary, we addressed thermal transport in metal- 
lic nanostructures contacted to the bulk of a 2DEG in 
quantizing magnetic fields showing that, thanks to the 
appearance of Landau gaps, a large heat current can 
flow off the N region. This leads to a significant elec- 
tron refrigeration. The performance offered over a wide 
range of temperatures makes this principle attractive as 
a solid-state quantum refrigerator in nanoscale metal- 
semiconductor systems. 
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